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We investigate the depinning transition for driven interfaces in the random-field Ising model 
for various dimensions. We consider the order parameter as a function of the control parameter 
(driving field) and examine the effect of thermal fluctuations. Although thermal fluctuations drive 
the system away from criticality the order parameter obeys a certain scaling law for sufficiently 
low temperatures and the corresponding exponents are determined. Our results suggest that the 
so-called upper critical dimension of the depinning transition is five and that the systems belongs 
to the universality class of the quenched Edward- Wilkinson equation. 

PACS numbers: 68.35.Rh,75.10.Hk,75.40.Mg 



I. INTRODUCTION 

Driven Interfaces in quenched disordered systems dis- 
play with increasing driving force a transition from a 
pinned interface to a moving interface (see e.g. Q and 
references therein). This so-called depinning transition 
is caused by a competition of the driving force and the 
quenched disorder. The first one tends to move the inter- 
face whereas the latter one hinders the movement. De- 
pinning transitions are observed in a large variety of phys- 
ical problems, such as fluid invasion in porous materials 
(see, for instance, Sec. 6.2 in 0] and references therein), 
depinning of charge density waves j^, 0], impurity pin- 
ning of flux-line in type-II superconductors j^], contact 
lines fl as well as in field driven ferromagnets, where the 
interface separates regions of opposite magnetizations . 

A well established model to investigate the depin- 
ning transition in disordered ferromagnets is the driven 
random-field Ising model (RFIM) (see for instance |t], [|, 
|9|, [If], |ll|). Here, the disorder induces some effective 
energy barriers which suppress the interface motion. A 
magnetic driving field H reduces these energy barriers 
but they vanish only if the driving field exceeds the crit- 
ical value H c . The transition from the pinned to the 
moving interface can be described as a continuous phase 
transition and its velocity v is interpreted as the order pa- 
rameter. Without thermal fluctuations (T = 0) the field 
dependence of the velocity obeys the power-law behavior 



v(h,T = 0) ~ hP 



(1) 



for h > 0, where h denotes the reduced driving field 
h = H/H c - 1. 
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The depinning transition is destroyed in the presence of 
thermal fluctuations (T > 0) which may provide the en- 
ergy needed to overcome local energy barriers. Although 
thermal fluctuations drive the system away from critical- 
ity the order parameter obeys certain scaling laws and for 
sufficiently low temperatures the order parameter can be 
described as a generalized homogenous function |12[ 



v(h,T) = A v(A~ 1/,/3 /i, \~^T), 



(2) 



similar to usual equilibrium second order phase transi- 
tions. Setting A~ 1///3 = 1 one recovers Eq. (Q) for zero 
temperature. Choosing A - ^ T — 1 one gets for the inter- 
face velocity at the critical field H c 



v(h = 0,T) ~ T 1 ^ . 



(3) 
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This power-law behavior was observed in two- and three- 
dimensional simulations of the driven RFIM [121 [l3| as 
well as in charge density waves in computer simulations 
and mean- field calculations [Q, [l4| . 

Furthermore, thermal fluctuations cause a creep mo- 
tion of the interface for small driving fields (H ^ H c ) 
characterized by an Arrhcnius like behavior of the veloc- 
ity [pl| |l6f . Recently, this creep motion was observed in 
experiments considering magnetic domain wall motion in 
thin films composed of Co and Pt layers fl7| , in renormal- 
ization group calculations |L8|, |l9| regarding the Edwards- 
Wilkinson equation with quenched disorder (QEW), as 
well as in numerical simulations of the RFIM |2(| . 

In equilibrium physics a scaling ansatz according 
to Eq. usually describes the order parameter as a func- 
tion of its control parameter and of its conjugated field. 
Although Eq. (^) can be applied to the depinning transi- 
tion, i.e. the temperature is a relevant scaling field, T is 
not conjugated to the order parameter. The conjugated 
field would support the interface motion independent of 
its strength. But strong thermal fluctuations destroy the 
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interface instead to support the interface motion. There- 
fore, one has to interpret the value of the thermal expo- 
nent if; carefully. For instance, it is not clear wether the 
obtained values of ip are a characteristic feature of the 
whole universality class of the depinning transition or 
just a characteristic feature of the particular considered 
RFIM. This point could be important for the interpreta- 
tion of experiments which naturally take place at finite 
temperatures. 

In this paper we reinvestigate the interface dynamics 
of the driven RFIM and focus our attention to higher 
dimensions d > 3. In particular we consider the scaling 
behavior at the critical point and determine the expo- 
nents (3 and ip. Our results suggest that the so-called 
upper critical dimension of the depinning transition of 
the driven RFIM is d c — 5. Above this dimension the 
scaling behavior is characterized by the mean-field ex- 
ponents. We compare our results with those obtained 
from a renormalization group approach of the quenched 
Edward- Wilkinson equation which is expected to be in 
the same universality class as the driven RFIM. A sum- 
mary is given at the end. 

II. MODEL AND SIMULATIONS 

We consider the depinning transition RFIM on cubic 
lattices of linear size L in higher dimensions (d > 3). The 
Hamiltonian of the RFIM is given by 

n = -\ Y. s - S i- H E 5 *-E^, (4) 

where the first term characterizes the exchange interac- 
tion of neighboring spins (Si = ±1). The sum is taken 
over all pairs of neighbored spins. The spins are coupled 
to a homogenous driving field H as well as to a quenched 
random-field h{ with (hi) — and (hihj) oc &y. The ran- 
dom field is assumed to be uniformly distributed, i.e., the 
probability p that the random field at site i takes some 
value hi is given by 

P[tlt) ~\ otherwise. (b) 

Using antiperiodic boundary conditions an interface is in- 
duced into the system which can be driven by the field H 
(see [jl2j for details). A Glauber dynamics with random 
sequential update and heat-bath transition probabilities 
is applied to simulate the interface motion (see for in- 
stance PH). 

Due to this algorithm not only spins adjacent to the 
interface but throughout the whole system can flip with 
a finite probability at temperatures T > 0. In general, 
this may cause nucleation which always starts with an 
isolated spin-flip. The minimum energy required for an 
isolated spin-flip is AE = 2(z J - H - A) with z = 2 d 
nearest neighbors on a bec lattice. Although the corre- 
sponding spin-flip probability l/[l + exp(A_E/T)] is small 



for the considered temperatures (T < 0.2), isolated spin- 
flips are possible and occur. But we observed that these 
spin- flips are unstable in our simulations, i.e., an isolated 
spin will flip back in the next update step. Thus the orig- 
inally induced interfaces is stable during the simulations. 

The moving interface corresponds to a magnetiza- 
tion M which increases in time t (given in Monte Carlo 
step per spin). The interface velocity, which is the basic 
quantity in our investigations, is obtained from the time 
dependence of the magnetization v = (dM/dt) where 
(. . .) denotes an appropriate disorder average. Starting 
with a flat interface we performed a sufficiently number of 
updates until the system reaches after a transient regime 
the steady state which is characterized by a constant av- 
erage interface velocity. 

As pointed out in previous works ^| an appropri- 
ate choice of the interface orientation is needed in order 
to recover that the interface moves for arbitrarily small 
driving field in the absence of disorder. An appropriate 
choice is to consider the interface motion along the di- 
agonal direction of a simple cubic lattice. For d > 3 an 
alternative is to examine the interface motion along the 
z-axis on a body-centered cubic (bec) lattice. Since it 
is much more convenient to implement the latter case in 
higher dimensions we consider in this work bec lattices, 
the more as the lattice structure usually does not affect 
the universal scaling behavior. 



III. D=3 

In the case of the three dimensional RFIM we consider 
bec lattices of linear size L < 250. A snapshot of a mov- 
ing interface in the steady state is presented in Fig.[l]. 
The obtained values of the interface velocities for T = 




FIG. 1: Snapshot of a moving interface at T = for L = 128, 
A = 1.7, and H = 1.37 in d = 3. In order to show the details 
we stretched the interface in the vertical direction by a factor 
15. 
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FIG. 2: Dependence of the interface velocity v on the driving 
field H for a bcc and simple cubic (sc) lattice, respectively. 
The inset shows v as a function of the reduced driving field h. 
The dash dotted lines are fits according to Eq. (w). 



are plotted in Fig. g. As one can see v tends to zero in the 
vicinity of H ss 1.36. Assuming that the scaling behav- 
ior of the interface motion is given by Eq. one varies 
H c until one gets a straight line in a log- log plot. Con- 
vincing results are obtained for H c = 1.357 ± 0.001 and 
the corresponding curve is shown in the inset of Fig.|^. 
For lower and greater values of H c we observe significant 
curvatures in the log- log plot (not shown). In this way 
we estimate the error-bars in the determination of the 



critical field. A regression analysis yields the value of the 
order parameter exponent (3 — 0.653 ± 0.026. This value 
agrees with j3 — 0.66 ± 0.04 which was obtained from a 
similar investigation |fl2"f where the interface moves along 
the diagonal direction of a simple cubic lattice (see inset 
of Fig.g). Furthermore, both results are in agreement 
with (3 = 0.60 ± 0.11 (h], [ll[], where in d = 3 the influ- 
ence of helicoidal boundary conditions in one direction 
and periodic ones in the other direction parallel to the 
interface was investigated on a simple cubic lattice. 

In order to determine the exponent ip we simulated 
the RFIM around the critical field for T = 0.025 n with 
n 6 {1, 2, 3, 4, 6, 8}. The obtained curves are shown in the 
inset of Fig. || According to Eq. (||) the interface velocity 
scales as 



(6) 



Plotting v{h,T)T- 1 ^ as a function of hT- 1 /^ one 
varies (3, ip, and H c until one gets a data collapse of the 
different curves. Convincing data collapses are observed 
for = 0.63 ± 0.06 tp = 2.33 ± 0.2 and H c = 1.360 ± 0.01 
and the corresponding curves are shown in Fig. ^. The 
obtained values of the order parameter exponent and of 
the critical field agree within the error-bars with the val- 
ues of the T = analysis. Furthermore our results are in 
agreement with similar investigations on a simple cubic 
lattice (/? = 0.63 ± 0.07 and ip = 2.38 ± 0.2, see ||). 



IV. D=4 

In order to determine the order parameter exponent 
of the four dimensional driven RFIM we measured the 
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FIG. 3: Scaling plot of the interface velocity for d — 3. The 
data are rescaled according to Eq. (^|). The inset shows the 
unsealed velocities for T = 0.025 n with n g {1,2,3,4,6,81 
(solid lines) in comparision to the T — data from Fig. H 
(dashed line). 



FIG. 4: The interface velocity of the four-dimensional model 
at T = 0. For sufficiently small fields the data obey a power 
law according to Eq. (|lj) (dotted dashed line) . For the fit we 
use only those data marked by filled symbols and we find 
H c = 1.258 ± 0.05 and f3 = 0.8 ± 0.06. 




FIG. 5: Scaling plot of the interface velocity for d = 4 and d - 
6, respectively (see inset). The data are rescaled according to 
Eq. 



interface velocity for bcc lattices of linear sizes L < 140. 
The obtained data for T = are shown in a log-log plot 
m Fig.|| After a transient regime which displays a finite 
curvature we observe an asymptotic power-law behavior 
for sufficiently small driving field H. A regression analy- 
sis yields /3 = 0.8 ± 0.06 and H c = 1.258 ± 0.002. 

To determine the exponent ip we simulated the RFIM 
in the vicinity of the critical field for T — 0.025 n where 
again n 6 {1,2,3,4,6,8} was choosen. Similar to the 
three dimensional case one varies the exponents as well as 
the critical field until one observes a data collapse. Good 
results are obtained for (3 = 0.73 ± 0.13, tp = 1.72 ± 0.27, 
and H c = 1.256±0.015. The corresponding data collapse 
is shown in Fig.^. Again, the obtained values of (i and 
H c confirm the above presented analysis for T = 0. 



V. D=5 

In the case of the five-dimensional RFIM system sizes 
from L = 10 up to L = 30 are simulated. Analyzing the 
interface motion at T = we observe that the velocity- 
field dependence can not be described by a pure power- 
law, i.e., Eq. (Q) fails. In Fig. |we plotted the logarithmic 
derivation of the velocity-field dependence 



A. 



d lav 
d \nh 



(7) 



which can be interpreted as an effective exponent. If the 
asymptotic scaling behavior obeys Eq. ([I]) the logarithmic 
derivative tends to the value of (3 for H — > H c . But as 
can be seen from Fig. 6 no clear saturation takes place for 
d — 5 as it is observed for the three and four dimensional 
model. The lack of a clear saturation could be explained 



FIG. 6: The effective exponents /3 e ff as a function of In h 
for various dimensions. The figure shows that the five- 
dimensional exponent does not display a clear saturation as 
the exponent of the lower dimensions do. The data of the 
six-dimensional model indicate that the expected saturation 
value /3 c g = 1 is reached significantly faster as compared to 
d — 5 (see text). The values of the two-dimensional system are 
obtained from fll3l. 



by a too large large distance from the critical point, but 
another reason is possible too. 

Significant deviations from a pure power law behav- 
ior [Eq. (|l|)] occur for instance at the upper critical di- 
mensional d c where the scaling behavior is governed by 
the mean-field exponents modified by logarithmic cor- 
rections. The scaling behavior around d c is well under- 
stood within the rcnormalization group theory (see for 
instance [^2|, |2^, |^|). For d > d c the stable fix point of 
the corresponding renormalization equations is usually a 
trivial fix point with classical mean-field exponents. This 
trivial fix point is unstable for d < d c and a different sta- 
ble fix point exists with nonclassical exponents. These 
exponents can be estimated by an e-expansion, for in- 
stance. For d = d c both fix points are identical and 
marginally stable. In this case the asymptotic form of 
the thermodynamic functions is given by the mean-field 
power-law behavior modified by logarithmic corrections. 
Applying this approach to the depinning transition, the 
corresponding ansatz reads 



v(h,T= 0) - h^MF | m /j| 



(8) 



where B denotes an unknown correction exponent. It is 
worth to mention that in contrast to the values of the 
critical exponents below the upper critical dimension the 
above scaling behavior does not rely on approximation 
schemes like e- or 1 /n-expansions ■ Within the renor- 
malization group theory it is an exact result in the limit 
h — > (see e.g. ]22], 26[ and references therein for RG 
investigations and [|27, p8[ for measurements). 
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FIG. 7: The rescaled interface velocity as a function of the 
driving field for d — 5. In order to display the logarithmic 
corrections we plot (v/h) 1 ^ vs. — ln/i [see Eq. (^)]. The solid 
line corresponds to the expected asymptotic scaling behav- 
ior for h — > [corresponding to In ft, — > ( — oo)] according to 
Eq.®. 



The value /3mf = 1 is reported for depinning transi- 
tions |2^, |30| . Thus we analyze v{h)/h as a function 
of | lnh\ and note again that Eq. (||) describes only the 
leading order of the scaling behavior, i.e., we expect that 
the asymptotic behavior of interface velocity obeys 

[v{h)/h} 1/B = const |ln/i|. (9) 

Therefore, we varied in our analysis the logarithmic cor- 
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FIG. 8: The scaling plot of the interface velocity v for the five 
dimensional model, i.e. at the upper critical dimension. The 
data are rescaled according to Eqs. ( |lo| , [ll"l ) using ip = 1.49, 
13 = 1, a = 0.22, and r = 0.19. 



rection exponent B and the critical field H c until we get 
this expected asymptotic behavior. The best results are 
obtained for B = 0.40 ± 0.09 and H c = 1.14235 ± 0.001 
and the corresponding scaling plot is shown in Fig.|^. 
The observed asymptotic agreement with Eq. (^|) corre- 
sponds to a logarithmically (l/|ln/i|) convergence of (3 c g 
to Pmf — 1) which explains why no clear saturation of 
the effective exponent could be observed for h — > in the 
five dimensional model. 

Similar to the T — scaling behavior one has to modify 
for T > the scaling ansatz since no data collapse could 
be obtained by plotting the data according to Eq. (g). 
Motivated by recently performed investigations of the 
scaling behavior of an absorbing phase transition around 
the upper critical dimension [ pl[ we assume that the scal- 
ing behavior of the order parameter obeys in leading or- 
der 



v(h,T) = T^mf | m T| CT v(x,l) 



(10) 



where the scaling argument x is given in leading order by 
x = hT-Ww^MF |lnT| r (11) 

with /?mf = 1< hi ° ur analysis we use the value V'mf 
I .-",'.) obtained from the analysis of the six dimensional 
RFIM (see next section). Therefore, we have to vary the 
exponents a and r in order to observe a data collapse ac- 
cording to Eqs. ( jl^ , |ll"| ). Convincing results are obtained 
for a = 0.22 ± 0.16, r = 0.19 ± 0.12 The corresponding 
data collapse is shown in Fig.||. 

VI. D=6 

Above the upper critical dimension the scaling behav- 
ior is characterized by the mean- field exponents, i.e., in 
leading order the interface velocity is given by 



v(h,T = 0) 



(12) 



In Fig. ^ we plot the velocity as a function of the driving 
field H obtained from simulations of system sizes L < 14. 
As one can see the velocity does not display the expected 
linear behavior. It seems that a linear behavior is only 
given for small velocities, i.e., our data do not display the 
pure asymptotic behavior [Eq. (|l2|)] . This is confirmed 
by the behavior of the effective exponent [Eq. ([?])] which 
increases fast for h — > but the actual saturation to /3 = 
1 does not take place for the considered values of h (see 
Fig. ||). To observe the asymptotic behavior one has to 
perform simulations closer to the critical point H c which 
requires to simulate larger system sizes. Unfortunately, 
the limited CPU power makes this impossible. 

An alternative is to take the curvature of the function 
v(h) into consideration and to assume that the leading 
corrections to the asymptotic behavior are of the form 



v(h,T = 0) = Vl h + v 2 h 2 + 0(h 3 



(13) 
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FIG. 9: The interface velocity as a function of the driving field 
for the six-dimensional model. The dashed line corresponds 
to a fit according to Eq. (|l3|). The inset displays v(h)/h as 
a function of the reduced driving field h. The resulting lin- 
ear behavior confirms that the deviations from the mean-field 
behavior can be described by quadratic corrections. 



which recovers Eq. ( |12| ) for h — ► 0. Fitting our data to 
this ansatz we get H c = 1.1537 ± 0.003, u x = 0.2546, 
and V2 = —0.0895. The corresponding curve fits the 
simulation data quite well as one can see from Fig.|^. 
In the inset of Fig.^j we plotted v(h)/h as a function 
of the reduced driving field h. According to the above 
ansatz [Eq. (13)] one gets a linear behavior, i.e., the de- 
viations from the pure mean-field behavior [Eq. (fl2"|) l can 
really be described as quadratic corrections. Thus we get 
that our numerical data are consistent with the assump- 
tion that the six dimensional RFIM depinning transition 
is characterized by the mean-field exponent /3mf = 1- 

Again we consider how thermal fluctuations affect 
the scaling behavior and analyse interface velocities ob- 
tained at different temperatures T = 0.025 n with n G 
{1, 2, 3, 4, 6, 8}. Similar to the situation below the upper 
critical dimension we assume that the scaling behavior 
of the interface velocity is given by Eq. (0) where the 
exponents are given by mean-field values. A convinc- 
ing data collapse is obtained for ipyiF — 1.49 ± 0.15 and 
H c = 1.153 ± 0.02 and is plotted in the inset of Fig. § 



VII. DISCUSSION 

A well established realization of interface pinning in 
a disordered media is the so-called quenched Edwards- 
Wilkinson (QEW) equation of motion which was inten- 
sively investigated in the last decade 0, |l8[ [l9| [2{| ||(], 
|32|. It is argued that the QEW equation as well as the 
driven RFIM are characterized by the same critical ex- 
ponents, i.e., both models belong to the same univer- 



FIG. 10: The critical exponents j3 and ij) as a function of 
the dimension. The solid line corresponds to the values of an 
e-expansion [Eq. ([14])]. 



sality class ]7|, [10|- Renormalization group analyses of 
the quenched QEW equation j2^, 30 predict, in accor- 
dance with d c — 5 and allow to estimate the critical 
exponents using an e-expansion. 



A recently performed 
two-loop renormalization approach yields [|30| 



Pqew 



0.040123 e 2 ±0(e 3 ) (14) 



where e denotes the distance from the upper critical di- 
mension, i.e. e = 5 — d (unfortunately, no error-bars can 
be estimated from an e-expansion). The corresponding 
values of the exponents as a function of the dimension 
are plotted in Fig. [l^. The numerically determined ex- 
ponents (3 of the driven RFIM (listed in Table |) are in 
a fair agreement with the values of the e-expansion. 

For the QEW equation the temperature exponent ip 
is not known. Therefore, a direct comparison with the 
obtained values of the driven RFIM is not possible. 



TABLE I: The exponents (obtained from simulations at 
T = and T > 0, respectively) and ij) of the depinning tran- 
sition of the RFIM for different dimensions. The values of 
the two dimensional model are obtained from Jl3[] . The criti- 
cal behavior at the upper critical dimension d c is additionally 
affected by logarithmic corrections. 

d /3 T =o Pt>o i> 



0.35 ±0.04 
0.653 ± 0.026 
0.80 ± 0.06 
1 
1 



0.33 ±0.02 
0.63 ±0.06 
0.73 ±0.13 
1 
1 



5.00 ±0.3 
2.33 ± 0.2 
1.72 ± 0.27 
1.49 

1.49 ± 0.15 
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VIII. CONCLUSIONS 

We studied numerically a field driven interface in the 
RFIM and determined the order parameter exponent j3 
as well as the temperature exponent ip. Below the up- 
per critical dimension d c — 5 the critical exponents de- 
pend on the dimension and the values of the exponents 
correspond to those of a two-loop renormalization group 
approach of the Edwards- Wilkinson equation j^o). This 
suggest that the depinning transition of the RFIM model 
belongs to the universality class of the quenched Edward- 
Wilkinson equation. At the upper critical dimension 
d c = 5 the scaling behavior is affected by logarithmical 
corrections. Above the upper critical dimension we ob- 



serve that the scaling behavior is characterized in leading 
order by the corresponding mean-field exponents. 
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